Optical frequency comb generation from a monolithic microresonator 
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Optical frequency combs[l], 0, H| provide equidistant frequency markers in the infrared, visible 
and ultra-violet 0, [a] and can link an unknown optical frequency to a radio or microwave frequency 
reference^, 0|. Since their inception frequency combs have triggered major advances in optical 
frequency metrology and precision measurements [H, 01 and in applications such as broadband 
laser-based gas sensingiSJ] and molecular fingerprinting 9j . Early work generated frequency combs 
by intra-cavity phase modulation^, [nj, while to date frequency combs are generated utilizing the 
comb-like mode structure of mode-locked lasers, whose repetition rate and carrier envelope phase 
can be stabilized f"P3]. Here, we report an entirely novel approach in which equally spaced frequency 
markers are generated from a continuous wave (CW) pump laser of a known frequency interacting 
with the modes of a monolithic high-Q microresonator 13] via the Kerr nonlinearity[3, EH- The 
intrinsically broadband nature of parametric gain enables the generation of discrete comb modes 
over a 500 nm wide span (~ 70 THz) around 1550 nm without relying on any external spectral 
broadening. Optical-heterodyne-based measurements reveal that cascaded parametric interactions 
give rise to an optical frequency comb, overcoming passive cavity dispersion. The uniformity of the 
mode spacing has been verified to within a relative experimental precision of 7.3 x 10 -18 . In contrast 
to femtosecond mode-locked lasers) 161] the present work represents an enabling step towards a 
monolithic optical frequency comb generator allowing significant reduction in size, cost and power 
consumption. Moreover, the present approach can operate at previously unattainable repetition 
rates [l7J exceeding 100 GHz which are useful in applications where the access to individual comb 
modes is required, such as optical waveform synthesis (T^i high capacity telecommunications or 
astrophysical spectrometer calibration[19l]. 



Optical microcavitiesjl^] are owing to their long tem- 
poral and small spatial light confinement ideally suited 
for nonlinear frequency conversion, which has led to a 
dramatic improvement in the threshold of nonlinear op- 
tical light conversion 2l( . In contrast to stimulated gain, 
parametric frequency conversion [22} does not involve cou- 
pling to a dissipative reservoir, is broadband as it does 
not rely on atomic or molecular resonances and consti- 
tutes a phase sensitive amplification process, making it 
uniquely suited for tunable frequency conversion. In the 
case of a material with inversion symmetry - such as silica 
- the non linear optical effects are dominated by the third 
order non linearity. The process is based on four-wave 
mixing among two pump photons (frequency vp) with a 
signal (ys) and idler photon (i/j) and results in the emer- 
gence of (phase coherent) signal and idler sidebands from 
the vacuum fluctuations at the expense of the pump field 
(cf. Fig.l). The observation of parametric interactions 
requires two conditions to be satisfied. First momentum 
conservation has to be obeyed. This is intrinsically the 
case in a whispering gallery type microcavity[20] since 
the optical modes are angular momentum eigenstates and 
have (discrete) propagation constants /3 m = ^ resulting 
from the periodic boundary condition, where the integer 
m designates the mode number and R denotes the cav- 
ity radius. Hence the conversion of two pump photons 
(propagation constant /3jv) into adjacent signal and idler 
modes (/3n-an, 0N+AN, AN = 1,2,3...) conserves mo- 



mentum intrinsically [14j (analogous reasoning applies in 
the case where the two annihilated photons are in differ- 
ent modes, i.e. for four- wave mixing, cf. Fig. lb). The 
second condition that has to be met is energy conser- 
vation. As the parametric process creates symmetrical 
sidebands with respect to the pump frequency (obeying 
hvi + hvs = 2/ii/p, where h is the Planck constant) it 
places stringent conditions on the cavity dispersion that 
can be tolerated since it requires a triply resonant cavity. 
This is a priori not expected to be satisfied, since the 
distance between adjacent modes ^fsr = \v m — ^m+il 
(the free spectral range, FSR) can vary due to both ma- 
terial and intrinsic cavity dispersion which impact n e s 
and thereby render optical modes (having frequencies 
where c is the speed of light in vacuo) 
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non-equidistant. Indeed, it has only recently been possi- 
ble to observe these processes in microcavities (made of 
crystalline [H CaF 2 and silica H). 
Importantly, this mechanism could also be employed to 
generate optical frequency combs: the initially generated 
signal and idler sidebands can interact among each other 
and produce higher order sidebands (cf. Fig. 1) by non- 
degenerate four- wave mixing (FWM)[24J which ensures 
that the frequency difference of pump and first order side- 
bands Av = \vp — vg\ = \vp — vi\ is exactly transferred 
to all higher order sidebands. This leads to an equidistant 
spectrum throughout the entire comb. This can be read- 
ily seen by noting that e.g. the 2 nd order sidebands are 
generated by mixing among the pump and first order sig- 
nal/idler sidebands (e.g. vgi = vp + vg — vi = vp — 2Av), 
which rigidly determines the spacing of any successively 
higher sidebands. Thus, provided the cavity exhibits low 
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FIG. 1: Broadband parametric frequency conversion from a 
75-/xm-diameter monolithic microresonator. Panel A: Spec- 
trum of the parametric frequency conversion observed in an 
75-/xm-diameter monolithic toroid microcavity when pumped 
with 60 mW continuous wave (CW) laser power at 1555 nm. 
The combination of parametric interactions and four-wave- 
mixing (FWM) gives rise to a broadband emission, spaced 
by the cavity free spectral range. Inset: Broadband para- 
metric conversion of a different sample generating more than 
70 parametric modes extending over a wavelength span of 
nearly 500 nm (launched power 130 mW). The asymmetry 
in the spectrum (with higher power in the higher wavelength 
sidebands) and the amplitude modulation of the emission is 
attributed to Raman amplification and variation of the ta- 
per fiber output coupling, respectively. Panel (B): Schematic 
of the processes that contribute to the parametric conversion: 
degenerate (left) and non-degenerate (right) four-wave-mixing 
among cavity eigenmodes. Momentum conservation is intrin- 
sically satisfied for the designated angular mode number (m) 
combinations. Panel (C): Scanning electron microscope image 
of a toroid microcavity on a silicon chip. 



dispersion, the successive four-wave mixing to higher or- 
ders would intrinsically lead to the generation of phase 
coherent sidebands with equal spacing, i.e. an optical 
frequency comb. Here, we report that microresonators 
allow realization of this process and generation of optical 
frequency combs. 

We employ ultra-high-Q monolithic microresonators in 
the form of silica toroidal microcavities[l3| on a silicon 
chip, which possess giant photon storage times (t) i.e. 
ultra-high quality factors (Q = 2ttvt > 10 s ) and small 
mode volumes. Highly efficient coupling is achieved us- 
ing tapered optical fibers [25|. Owing to the high cir- 
culating power, parametric interactions are readily ob- 
served at a threshold of approx. 50 /iW. When pump- 
ing with a continuous wave (CW) 1550-nm laser source, 



we observe a massive cascade and multiplication of the 
parametric sidebands extending to both higher and lower 
frequencies. Fig. la shows a spectrum in which a 
75-/im-diameter microcavity was pumped with 60 mW 
power, giving rise to an intra-cavity intensity exceed- 
ing 1 GW/cm 2 . The parametric frequency conversion 
could extend over more than 490 nm (cf. Fig la inset), 
with the total conversion efficiency being 21.2 % (The 
highest observed conversion efficiency was 83 % by work- 
ing in the over-coupled regime [II]). These bright side- 
bands (termed Kerr combs in the remaining discu ssions) 
could be observed in many different samples. Also, in 
the largest fabricated samples (190 /im diameter) 380- 
nm broad Kerr combs comprising 134 modes spaced by 
375 GHz could be generated at the expense of slightly 
higher pump power (cf. Supplementary Information). 

To verify that the Kerr comb indeed contains equidis- 
tant frequencies, we employed a fiber-laser based optical 
frequency comb_26| from Menlo Systems (termed "refer- 
ence comb" in the remaining discussion) as a reference 
grid whose repetition rate is / rep = 100 MHz. The prin- 
ciple underlying our measurement is that the beating 
generated on a photodiode by superimposing the refer- 
ence comb with the Kerr comb will produce beat notes 
which constitute a unique replica of the optical spectrum 
in the radio frequency domain, provided that the highest 
produced beat note in the detection process is < / re p/2 
(cf. Fig. 2b), similar to multi-heterodyne frequency 
comb spectroscopy (27). Specifically, if the Kerr comb is 
equidistant, the beat notes with the reference comb will 
constitute an equidistant comb in the RF domain (with 
frequency spacing A/, where A/ = (Av mod / re p))- 
Fig. 2a shows the experimental setup for the optical 
beat measurement. In brief, an external cavity laser at 
1550 nm was used as pump laser (cf. Fig 2 main panel). 
The Kerr lines of the microcavity were superimposed 
with the reference comb in a beat note detection unit 
(BDU), consisting of polarizing optics to combine the 
reference and Kerr combs and a grating to select the 
desired region of spectral overlap. In this manner, the 
beating of 9 simultaneously oscillating parametric modes 
(covering > 50 nm of wavelength span) were recorded, 
as shown in Fig. 2c. Remarkably, from the equidistant 
spacing of the radio-frequencies, it is found that the 
generated sidebands are equidistant to within less than 
5 kHz (as determined by the measurement time of 200 

(IB). 

To improve the accuracy, we developed an additional 
experiment where we measured the beat-notes of three 
Kerr modes with the fiber-reference comb using three 
separate BDUs (cf. Fig. 3a) each counting a single beat 
of radio frequency (/o, /i, /a)- A signal-to- noise ratio 
exceeding 30 dB in 500 kHz bandwidth was achieved, 
sufficient to use radio-frequency counters, which were 
all referenced to a 10 MHz reference signal provided by 
the MPQ hydrogen Maser. The beat-note measured on 
BDU1 was used to implement an offset lock between a 
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FIG. 2: Parametric beat note setup. Panel (A) shows the ex- 
perimental setup consisting of an external cavity laser (ECL) 
coupled to an ultra-high-Q monolithic microresonator in a 
nitrogen environment via a tapered fiber. The parametric 
output is coupled into one port of a beat note detection unit 
(BDU). The second port of the BDU is coupled to a mode- 
locked femtosecond erbium doped fiber laser that serves as 
a reference comb. A grating is used to select a spectral re- 
gion of the Kerr modes and a PIN Si photodiode records their 
beatings with the reference comb (See supplementary infor- 
mation for details). Panel (B) shows the measurement prin- 
ciple. The beating of the reference comb with the parametric 
lines yields beat frequencies in the radio-frequency domain. 
Panel (C): RF Spectrum of 9-simultaneously oscillating Kerr 
modes, exhibiting a uniform spacing. 



single reference comb line and the pumping laser by a 
known offset frequency (fo)- The second (third) counter 
measured the N th (M th ) mode of the Kerr comb as 
shown in Fig. 3. For equidistant mode spacing, the 
second (third) BDU gives rise to the beat frequency 
h = fo + A/ X N (f 2 = f + Af X M) . The uniformity of 
the Kerr comb was then checked by comparing the varia- 
tion in the mode spacing, i.e. e = ~ ^a/" ■ Alterna- 
tively, direct counting of the ratio (/i — fo) / (/ 2 — fo) — 
jj was implemented (using frequeny mixing and ratio 
counting, cf. SI). Fig. 4b shows the result of this mea- 
surement (for N = 5, M = 7) and a counter gate time 
(r) of 1 second and more than 3000 records. The scatter 
in the data follows a Gaussian distribution (and follows a 
dependence of the Allan deviation, cf. Fig. 4b). The 
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FIG. 3: Probing the equidistance of the comb structure. Panel 
(A): Simplified schematic of the setup which consists of three 
beat note detection units (BDUs) to measure the beating of 
three Kerr modes with the fiber based reference comb. All 
BDUs are referenced to the MPQ hydrogen Maser as fre- 
quency standard. The first BDU is used to implement a phase 
lock between one comb line of the reference comb and the 
pump laser (which constitutes one mode of the Kerr comb). 
Panel (B): The parametric spectrum under consideration for 
validating the equidistance of the comb modes. 



cavity modes of this measurement span over approx. 21 
nm and yield a deviation from the mean of ±5 mHz. Note 
that the wavelength span that could be used for the mea- 
surement is currently limited by the gain bandwidth of 
an EDFA, which had to be used to amplify the reference 
comb to have sufficient power to run three BDUs simul- 
taneously. The results for different gate times and for the 
two different counting methods are shown in Table 1 (the 
complete list is contained in the supplementary informa- 
tion). The weighted average of these results verifies the 
uniformity of the comb spacing to a level of 7.3 x 10~ 18 
(when referenced to the optical carrier). Normalized to 
the bandwidth of the measured Kerr lines (2.1 THz), this 
corresponds to 5.2 x 10~ 16 . This accuracy is on par with 
measurements for fiber based frequency combs [26j and 
confirms that the generated Kerr combs exhibit uniform 
mode spacing. 



Next we investigated the role of dispersion under- 
lying the observed comb generation. Dispersion in 
whispering-gallery-mode (WGM) microcavities is char- 
acterized by the deviation in the free spectral range 

A^FSR = (i-m+l — V m ) — (v m — V m -\) = h> m+ i — V m -\ — 

2v m and has two contributions. Geometrical disper- 
sion accounts for a negative FSR dispersion, given by 



AfFSR ~ -0.4U 



m 



_5 / 3 where R the cavity radius 



'2-K-n-R 

(cf. supplementary information). Material dispersion on 
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TABLE I: Summary of the experimental results on the accu- 
racy of the mode spacing. The weighted mean of all measure- 
ments (including SI) yields a relative accuracy of 7.3 x 10~ 18 . 



250 
200 

! 150 

i 

; 100 

50 



Total Counts = 3493 

Mean = -0.91 mHz ± 5.46 mHz 

Standard Deviation = 322 mHz, 



i 




Gatetime = 1 s 



-1.0 -0.5 0.0 0.5 1.0 

Variation of the Mode Spacing (Hz) 



= 1E-15 



£ 1E-17 









* Measured Allan Deviation 








Fit 








<r =1.29e-15T'°' 53 

A 















10 100 
Gatetime t (s) 



FIG. 4: Frequency counting experiment. Panel (A): The de- 
viation from equidistant mode spacing for a gate time of 1 
second for the parametric spectrum and measurement setup 
described in Fig 3. For this measurement 3493 points were 
recorded. The solid red line is a Gaussian fit to the distri- 
bution. The deviation from the mean implies an accuracy 
of the mode spacing at the mHz level, confirming the comb 
like structure of the parametric spectrum. Panel (B): Allan 
deviation as a function of gate time, exhibiting an inverse 
square-root dependence on the gate time, as determined by 
the fit (solid red line). 
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the other hand is given by A^fsr 
where GVD = — is the group velocity dispersion 
parameter. Since the GVD of silica is positive for wave- 
length greater than 1.3 /xm (anomalous dispersion), it 
can compensate the intrinsic resonator dispersion (caus- 
ing AfpsR > 0). Indeed we measured a positive disper- 
sion (cf. SI) which equates to only ca. 20 MHz over a 
span of ca. 60 nm. This low value indicates that the 
present experiments are carried out close to the zero dis- 
persion wavelength, in agreement with theoretical pre- 
dictions. 

Note that the residual cavity dispersion exceeding the 
"cold" cavity linewidth does not preclude the paramet- 
ric comb generation process. This can be explained in 
terms of a nonlinear optical mode pulling process as re- 
ported in Ref. [l4| ■ The strong CW pump laser will in- 



duce both self-phase modulation (SPM) and cross-phase 
modulation (XPM)[28(, the latter being twice as large 
as the former. The resultant XPM and SPM induced 
refractive index changes will shift the cavity resonance 
frequencies by different amounts, thereby causing a net 
change in the (driven) cavity dispersion from its passive 
(un-driven) value (l4|. This nonlinear mode pulling can 
provide a mechanism to compensate the residual cavity 
dispersion. 

Regarding future experimental work in light of applica- 
tions in metrology, we note that absolute referencing can 
be attained by locking the pump laser to a known atomic 
transition and locking the mode spacing to a microwave 
reference (such as a Cs atomic clock). The latter requires 
that the two degrees of freedom of the comb, its repeti- 
tion rate (i.e. mode spacing, Av) and frequency offset, 
i.e. ^ceo = { v o mod Av) to be controlled independently. 
Indeed it could already been shown in a proof of concept 
experiment that it is possible to lock two modes of the 
Kerr comb simultaneously to two modes of the reference 
comb, showing that we are able to control both vceo and 
Av. The two actuators used for this lock are the detun- 
ing of the pump laser from the microcavity resonance and 
the pump power, which affects the optical pathlength of 
the cavity via the thermal effect and the nonlinear phase 
shift. 

Pertaining to the implications of our work, we note that 
the present observation of a monolithic frequency comb 
generator could potentially prove useful for frequency 
metrology, given however further improvements. Evi- 
dently a readily measurable repetition rate would prove 
useful when directly referencing the optical field to a mi- 
crowave signal|2[. To this end a 660-/mi-diameter micro- 
cavity would already allow operating at repetition rates 
< 100 GHz, which can be detected using fast photo- 
diodes. On the other hand, a large mode spacing as 
demonstrated here could prove useful in several appli- 
cations, such as line-by-line pulse shaping, calibration of 
astrophysical spectrometers or direct comb spectroscopy. 
The high repetition rate from an on chip device may also 
prove useful for the generation of multiple channels for 
high capacity telecommunications (spacing 160 GHz) and 
for the generation of low noise microwave signals. Fur- 
thermore, we note that parametric interactions do also 
occur in other types of microcavities - e.g. CaF2[15] - 
provided the material exhibits a third order nonlinearity 
and sufficiently long photon lifetimes. As such the cavity 
geometry is not conceptually central to the work and the 
reported phenomena should become equally observable 
in other types of high-Q microresonators, such as silicon, 
SOI or crystalline based WGM-resonators. Indeed the 
recent observation of net parametric gain[2^| on a silicon 
chip is a promising step in this direction. 
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I 

I. GENERATION OF KERR COMBS AT LOWER REPETITION RATES 



Figure O shows the Kerr comb spectrum at a lower repetition rate mentioned in the main paper. The repetition 
rate is 375 GHz, corresponding to a free spectral range of 3 nm. With larger samples it should be possible to generate 
repetition rates smaller than 100 GHz which permits the direct measurement of the repetition rate with high-bandwidth 
photodiodes. 
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FIG. 5: Kerr comb generated in an 177-/im-diameter toroid. The total power in the spectrum (pump line + generated sidebands) 
is around 500 mW distributed over more than 134 lines. The free spectral range is 3 nm. 



II. BEAT NOTE EXPERIMENTS BETWEEN THE FIBER LASER COMB AND KERR COMB 



To demonstrate the equidistant nature of the parametric Kerr lines, a reference frequency comb in the form of a 
mode locked erbium fiber laser is used (from Menlo Systems GmbH). The principle underlying the measurement is 
similar to the concept of multi- heterodyne spectroscopy [13] • Assuming that the reference comb produces a spectrum 
with frequencies /ceo + ?i • /rep (where / rep is the repetition rate, / cco is the carrier envelope offset frequency and n is an 
integer number of order 2 • 10 6 ) and the Kerr comb produces frequencies uq + m ■ (to integer), the signal generated 
by interfering the two combs will have an imprinted radio frequency (RF) beat note spectrum. If the reference comb's 
repetition rate is adjusted such that a multiple of it is close to the Kerr mode spacing, i.e. too 1 /rep ~ Av (with an 
integer Too), then the N different Kerr comb lines will generate N different RF beat notes which will again be evenly 
spaced, i.e. their frequencies are /o + A/ • k (with A/ = (Ai/mod/ rep ) and k integer). 

The experimental setup is depicted in Figure 2 of the main paper. A tunable external cavity diode laser (ECDL) 
is used to pump a microtoroid resonance as detailed in [30l ] and [lij]. Since the cavity resonances are polarization 
dependent, a in-fiber polarization controller is used to adjust the polarization of the pump laser. The microtoroid 
is placed in a sealed enclosure containing a nitrogen atmosphere, to avoid the deposition of water on the surface of 
the silica toroid which has strong absorption bands in the 1550-nm regime. In the microresonator a spectrum of 
modes is generated via nonlinear parametric interactions and four- wave mixing (see main paper). The output signal 
of the tapered optical fiber (containing the parametric modes that are outcoupled from the microresonator back to 
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the tapered fiber) is split by two 3 dB couplers and monitored with a photodiode connected to an oscilloscope and 
an optical spectrum analyzer. Another fraction of the taper output is sent to a "beat detection unit" (BDU) and 
superimposed with a fiber-laser-based reference frequency comb with a repetition rate of 100 MHz[2^|. The BDU 
consists of quarter wave plates and half wave plates to prepare orthogonal linear polarization in the two input beams, 
which are subsequently combined using a polarizing beam splitter. By means of a half-wave plate, an adjustable linear 
combination of the two input beams' polarizations is then rotated onto the transmission axis of a polarizer, where the 
two input beams interfere. To increase the signal-to- noise ratio (SNR), the spectral region containing the Kerr comb 
lines is selected by a grating and finally detected with a PIN InGaAs photodiode (Menlo Systems FPD 510). An 
oscilloscope with a built-in FFT routine is utilized to analyze the radio frequency spectrum. For rough analysis an 
electronic spectrum analyzer is used. Since the repetition rate of the reference comb is around 100 MHz the beat note 
frequencies between a laser line and the reference comb are in the range of MHz to 50 MHz. Now the repetition rate 
of the reference comb is adjusted until (A^mod/ rep ) is a small frequency such that for all k of interest the condition 
< /o + k ■ A/ < / r ep/2 is fulfilled. The observation of an equidistant RF beat "comb" then provides proof for the 
equidistance of the Kerr comb. 



III. MEASURING THE ACCURACY OF THE MODE SPACING USING HETERODYNE 

SPECTROSCOPY 

A. Measuring with two counters 

To verify the equidistance of the Kerr comb modes it is necessary to know the frequencies of three Kerr comb 
modes simultaneously. The frequency counting is achieved by using radio frequency counters that are connected to a 
photodiode in a beat note detection unit (cf. figure 2 in the main paper). To determine the frequencies of three Kerr 
comb modes at the same time, three beat note detection units (BDU) have been built. By tuning the grating of the 
BDUs it is possible to measure the beat note frequency of a single Kerr comb line with a reference comb line. For 
simplicity reasons, one BDU is used to lock the diode laser pumping the microcavity to a single mode of the reference 
comb. Additionally the repetition rate of the reference comb is locked to a frequency of around 100 MHz, stabilized 
by a 10 MHz frequency standard generated by an in-house hydrogen maser. The two remaining beat detection units 
are placed at the output of the microcavity and the gratings are adjusted in a way that each of them counts a different 
Kerr comb mode. Note that the output of the reference comb had to be amplified with an EDFA to obtain sufficient 
power to run three BDUs simultaneously (a single line of the reference comb contains ca. 10 nW optical power). With 
this setup it was possible to achieve signal-to-noise ratios for the Kerr sideband beat notes of more than 30 dB at a 
resolution bandwidth of 500 kHz (Additional RF filters with a 3-dB-bandwidth of 3 MHz have been used to filter out 
background noise). In the present experiment we focused on counting the 5 th (beat note frequency /i) and the 7 th 
Kerr comb sideband (beat note frequency /2X whereas the pump laser was phase locked to the fiber laser reference 
comb such that its beat with the reference comb was fixed to a frequency /q. Note that the pump laser already 
constitutes one tooth of the Kerr comb. For an equally spaced Kerr comb we therefore expect fi = fo + N- A/ and 
h = fo + M ■ A/ with N = 5 and M — 7 to be the beat note frequencies of the sidebands. The variation of the mode 
spacing e of the Kerr comb is given by 

H — fi fi — fo / 1 s 

e = , (1 

M — N N w 

which is zero for an equally spaced comb. With the measured values for f\ and fi and the known frequency fo it 
is possible to calculate the variation of the mode spacing e. The two counters are referenced to the same frequency 
standard as the offset lock for the pump laser and are externally triggered with a signal from a pulse generator. This 
external triggering was necessary since the mode spacing of the Kerr comb was fluctuating by approximately 40 kHz 
r.m.s., giving rise to "breathing" of the Kerr comb modes(cf. figure^. Hence, it proved cricital for a high accuracy 
that the two counters measured simultaneously, to allow the cancellation of the common fluctuations. 



B. Measuring the ratio of the distance to the sidebands 

To avoid the synchronization problems mentioned before, the experimental setup depicted in figure [7| was used. In 
brief, the three counter signals were first electronically mixed with f and filtered yielding only the distance between 
pump and the N th (M th ) Kerr-sidebands. With this setup, a slightly smaller standard deviation of the measurements 
could be achieved by using just one counter with two inputs to measure the ratio R of the distance between the pump 
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FIG. 6: Drift of the mode spacing of the Kerr comb when not stabilized. The mode spacing exhibits fluctuations of approximately 
40 kHz for short time scales and some slower thermal drifts for time scales of several minutes. Note that these fluctuations of 
the mode spacing do not affect the equidistance of the modes of the Kerr comb since they are fluctuating simultaneously. 



beat and the two sideband beats, 



R 



fa — fo 
fi - fo 



(2) 



Solving this for fa and using equation [T] we obtain the dependence of the variation of the mode spacing e from the 
ratio R: 



fx - fo D . , , r - ( 1 
e = — -R+(fo- h. 



M-N 



1 



M-N N 



(3) 



Using the frequency difference fi — fa, which was set to approximately 10 MHz, it is possible to derive the variation 
of the mode spacing e by measuring the frequency ratio R. 



IV. EXPERIMENTAL RESULTS OF THE COUNTER MEASUREMENTS 



Table |TT] shows the experimental results from the measurements of the Kerr comb equidistance. Note that a total of 9 
data points out of the 8382 measurements from tablcHTlhave been removed from analysis. These data points have been 
separated by the other data points of the respective measurements by more than 15 standard deviations. Assuming 
a Gaussian distribution (which was indeed found for the remaining 8373 measurements) the probability of measuring 
a point 15 standard deviations off as given by the cumulative error function is (1 - erf (15/^2)) « 7.3 x 1CT 51 . These 
points are believed to originate from some local perturbations in the lab leading to a temporary reduction of the 
signal-to-noise level of the radio frequency beat notes. The weighted mean e w in table |TT] has been calculated with the 
inverse squared standard error of the mean as weight: 



(4) 



The weighted mean calculated from all measurements leads to a variation of the modespacing of e w — —0.8 mHz ± 
1.4 mHz. Normalized to the optical carrier frequency of 192 THz, this leads to an accuracy of the equidistance of 
7.3 x 10~ 18 . 
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FIG. 7: Panel A). Experimental setup to measure the ratio of the frequency separation between pump laser and two different 
Kerr comb sidebands. ECDL = External Cavity Diode Laser. Beat note detection unit f (BDU1) is used to phase lock the 
pump laser line from the Kerr comb to the reference comb with an offset frequency /o. BDU2 (BDU3) is adjusted to measure 
the beat note frequency between the N th (M* h ) Kerr comb line and the reference comb. By mixing these frequencies down 
with the offset lock frequency /o using electronic mixers, new frequencies TV ■ Af and M ■ Af are generated. The ratio of these 
frequencies is M/N = 1.4 for the 7 th and the 5 th sideband. Panel B) shows a measurement of the frequency ratio of the radio 
frequency beat notes generated from the 7 th and the 5 th sideband of the Kerr comb. 



V. MEASUREMENT OF CAVITY DISPERSION 



To measure cavity dispersion, we employ the arrangement shown in figure[5J In brief, we first lock an external cavity 
laser around 1550 nm to one of the fundamental WGM cavity modes (the same resonance that gives rise to cascaded 
sidebands at higher power). The cavity resonance of the monolithic microresonator is locked to the external cavity 
laser by virtue of the thermal self locking technique 31]. The power is chosen to be far below the parametric threshold 
< 85 fiW but sufficient to entail a stable lock. Next, the frequency comb is offset- locked to the external cavity laser 
by recording the beat note signal in a separate beat note detection unit (for working principle of the beat detection 
unit see last section). To achieve stable locking the generated beating is filtered and amplified yielding a SNR of 
ca. 25 — 30 dB (at a resolution bandwidth of 400 kHz). For dispersion measurement the frequency comb must be 
locked at an arbitrary detuning with respect to the ECDL. The latter is accomplished by mixing the beat note with a 
(variable) reference signal (/offset) down to 10 MHz and implementing a phase lock with feedback on the fiber comb's 
repetition rate (/ re p) by controlling the cavity length using a mirror mounted on a piezoelectric tube (Note that all 
RF generators and analyzers are stabilized using an in- house 10-MHz- reference). Owing to the fact that the cavity 
linewidth is < 5 MHz and the repetition rate of the fiber comb (FC) is 100 MHz, not more than one FC comb mode 
at a time can be resonant with one microresonator mode. Since measuring the coupling of an individual comb mode 
into the resonator in transmission is difficult, we measure the reflection of the cavity induced by modal coupling 32]. 
By variation of /offset (and by recording simultaneously / re p) this allows to resolve the linewidth of individual cavity 
modes in reflection when using an OSA in zero-span mode. Hence this measurement provides an accurate means to 
measure frequency gap (free spectral range) between two cavity resonances v m and v m +Am modulo the repetition rate 
of the fiber comb {{v m — ^ m +Am)niod/ r ep). The low power of the individual FC lines (ca. 10 nW) ensures that the 
probed cavity mode is not thermally distorted. To remove the ambiguity in the number of comb lines (n) between the 
FSR of the cavity i.e. n = \_(v m — ^ m +Am)//repJ a second measurement was carried out with a different repetition 
rate, which allowed to retrieve n. So the actual free spectral range between two cavity resonances can be derived by: 



^FSR — /beatnote + Tl • / rC p 
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Gate time (s) 


Readings 


Mean Value for e (mHz) 


StdDev of e (Hz) 


Counting Method 


0.03 


217 


-33 ± 556 


8.2 


ratio 


0.1 


223 


-80 ± 181 


2.7 


ratio 


0.3 


293 


2.4 ± 50.1 


0.86 


ratio 


1 


3493 


-0.91 ± 5.46 


0.32 


2 counters 


1 


3499 


3.9 ± 10.1 


0.60 


2 counters 


1 


98 


-40.1 ± 27.4 


0.27 


ratio 


1 


179 


8.0 ± 25.5 


0.34 


ratio 


3 


173 


5.8 ± 12.6 


0.17 


ratio 


10 


22 


-17.9 ± 15.0 


0.070 


ratio 


30 


39 


1.65 ± 7.41 


0.046 


ratio 


60 


72 


-1.88 ± 3.00 


0.025 


ratio 


100 


18 


1.12 ± 5.98 


0.024 


ratio 


100 


42 


-0.26 ± 2.69 


0.017 


ratio 


300 


14 


-0.82 ± 2.83 


0.011 


ratio 


Weighted Mean e m : 


-0.8 mHz ± 1.4 mHz 







TABLE II: Complete list of the Kerr comb measurements with frequency counters. StdDev = standard deviation of the 
distribution. The last column shows the used method to acquire the data: "2 counters" stands for the measurements with two 
externally triggered counters (one for each Kerr sideband) and "ratio" stands for the method with one counter that directly 
measures the ratio between the distance between pump and two different Kerr comb lines (both methods are explained in the 
preceding section). As expected, the standard deviation of the measurements reduces with increasing gatetime. The total 
measurement time is 6 h 37 min. 



Figure [5] shows the experimental result of the dispersion measurement. The used cavity had a free spectral range 
(FSR) of 7.9 nm, which corresponds to 0.96 THz. Plotted in figure [H is the accumulated dispersion of the FSR, 
which we express for convenience as {v m +i — v m ) — {v\ — v$). Here, the v m are the resonance frequencies of a "cold" 
microcavity. For this measurement, i>q is a resonance at 1585 nm (189 THz). From the graph it can be derived that 
the accumulated dispersion is 2.6 MHz per FSR (i.e. positive dispersion). 



VI. DISPERSION PREDICTIONS 



The dispersion in our microcavities has two contributions. First, whispering-gallery mode microcavities exhibit 
an intrinsic variation of the free spectral range owing to the resonator geometry. The resonance frequency of the 
fundamental mode of a microsphere is approximately given by [33| 

,+ 1/2 + fft -i- +... , (6) 



2nnR 



where c is vacuum light speed, n the refractive index, R the cavity radius and r\\ the first zero of the Airy function 
sa 2.34). Hence, the variation of the free spectral range 



is given by 



d 2 v m 

A^FSR = {v m +l ~ V m ) — (v m — V m -l) = Vm+l + U m -i — 2v m w — — — (7) 



A , FSR = . 2L ( m + 1 / 2 ) ~ 5/3 „ _ .41^m- 5 /3 < o (8) 

tbK 2nnR 18 V 2 J 2%nR V ; 



Evidently, the free spectral range reduces with increasing frequency corresponding to a negative group velocity disper- 
sion (GVD), i. e. low frequency modes exhibit a shorter round trip time than high frequency modes. Supplementary 
figure S3 shows the variation for a 40- and 80-micron-radius microsphere. 

A second contribution comes from the dispersion of the fused silica material constituting the resonator. Its contri- 
bution can be estimated by considering that the refractive index n is actually a function of frequency (and therefore 
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FIG. 8: Experimental setup of the dispersion measurement. The beat detection unit on the lower left side is used to establish 
an offset lock between the external cavity diode laser (ECDL) and the fiber laser frequency comb. Therefore the signal from 
the photodiode in the beat detection unit is first filtered with a 50 MHz lowpass to remove the strong signal of the 100 MHz 
repetition rate of the fiber laser comb. Subsequently the beat note signal is mixed down to 10 MHz with a variable frequency 
generator (10. .60 MHz) and compared with a stable 10 MHz RF reference. The output of the comparator is sent to a PI 
feedback amplifier which is connected to a piezo-mechanical control of the repetition rate of the fiber laser. By adjusting the 
variable frequency generator one can change the distance between the laserline of the ECDL and the next comb line to an 
arbitrary value between MHz and / rep /2. The ECDL and the fiber comb are furthermore coupled to the microcavity with 
a microtoroid resonance thermally locked to the ECDL. To measure the distance between two cavity resonances an optical 
spectrum analyzer (OSA) in zero span mode is set to a wavelength of a different cavity resonance than the one pumped by the 
ECDL. Next, the offset lock is changed until a reflection signal of the fiber comb is detected on the OSA. Once this is achieved 
the ECDL and one mode of the fiber comb are on resonances with two different modes of the microcavity. This means the FSR 
can be derived as /beatnote + n ■ / rcp . 



mode number m), n = n(m). Neglecting geometric dispersion, the GVD of fused silica alone would lead to a FSR 
variation of 

A " FSR w (^Ri? ' m ) * 4^w ' GVD ' (9) 

where 

GVD = -^ (10) 



c 



zdX 2 



is the group-velocity dispersion of fused silica. This material parameter is well-known to change its sign in the 
1300-nm wavelength region from about — lOOps/nm • km at 800 nm to +20ps/nm • km at 1550 nm. Combining the 
two contributions, the positive sign of the GVD allows us in particular to cancel the geometric dispersion of our 
resonators to some extent, rendering the FSR nearly constant over a wide frequency span. Figure [TU] displays the FSR 
variation for an 80- and 160-micrometer diameter microsphere, considering both material and geometric dispersion. 
Importantly, a zero dispersion point close to our operating wavelength occurs. Note that for a toroidal microcavity the 
location of the zero dispersion point is expected to be shifted to shorter wavelengths owing to the different resonator 
geometry. This expectation is borne out of finite element simulations showing that the resonance wavelength for a 
given m value is shorter in a microtoroid cavity as compared to a microsphere [34j . 
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FIG. 9: Dispersion measurement of an 80-/im-diameter monolithic microresonator. The figure shows the accumulated variation 
(i.e. dispersion) of the free spectral range i.e. (v m +i — u m ) — (yi — vo). The variation of the FSR at higher frequencies (shorter 
wavelength) is referenced to the free spectral range recorded between 1577 nm (u\) and 1584 nm (vq). The shaded region 
denotes experimental uncertainty, the dotted line denotes a linear fit. As expected for a whispering- gallery mode dominated 
by material dispersion, the free spectral range increases for shorter wavelength. 
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FIG. 10: Variation of the free spectral range of a whispering-gallery microsphere resonator (i.e. Az^fsr = fm+i + v m -\ — 2v m )- 
Shown is the FSR dispersion for two resonator radii (40 /im and 80 fim) including the effect of silica dispersion via the Sellmeier 
equation. Resonance locations were calculated using an asymptotic expansion of the microsphere resonance locations. Owing 
to the different signs of silica material and resonator dispersion, a zero dispersion point exists in the infrared. 



APPENDIX A: SYMBOLS USED THROUGHOUT THIS WORK 
Symbols Designation 



v m Optical microcavity mode (with angular mode numbcrm) 

fpsR Optical microcavity free spectral range (z^fsr = \ v m ~ ^m+il) 

A^fsr Optical microcavity variation of the free spectral range (A^pgR = v m +i + v m -i ~ 

v cco Kerr comb carrier envelope offset frequency 

Av Kerr comb mode spacing 

/ rcp Fiber reference comb repetition rate 

/coo Fiber reference comb carrier envelope frequency 

/o,i,2 Beat note unit (BDU) frequencies 

A f Frequency spacing of the multi- heterodyne beat comb 
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